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Fast control of topological vortex formation in BEC by connter-diabatic driving 
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Topological vortex formation has been known as the simplest method for vortex formation in 
BEC of alkali atoms. This scheme requires inversion of the bias magnetic field along the axis 
of the condensate, which leads to atom loss when the bias field crosses zero. In this Letter, we 
propose a scheme with which the atom loss is greatly suppressed by adding counter-diabatic magnetic 
field. A naive counter-diabatic field violates the Maxwell equations and we need to introduce an 
approximation to make it physically feasible. The resulting field requires an extra currents, which 
is experimentally challenging. Finally we solve this problem by applying a gauge transformation so 
that the counter-diabatic field is generated by controlling the original trap field with the additional 
control of the bias field. 
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Demonstration of quantized vortices in a superfluid 
is a manifestation of the nonvanishing order parameter. 
Topological structure of vortices often reflects the mani¬ 
fold of the order parameter space. It is, therefore, natural 
to seek for the method to experimentally demonstrate 
formation of vortices once BEC of alkali metal atoms 
has been realized. Topological method of vortex forma¬ 
tion is the simplest one compared with other proposals, 
where the spinor structure of the order parameter is fully 
utilized [ll-l7|. Although the vortex [6| and monopole 
d, 1^ in BEC have been observed experimentally, estab¬ 
lishment of methods to create such interesting quantum 
states more accurately in time much shorter than coher¬ 
ence time is important for further experimental study 
and applications. 


A condensate in the weak held seeking state (WFSS) 
is trapped in a quadrupole magnetic held in the a;y-plane 
with a uniform bias held along the 2 ;-axis hrst and, sub¬ 
sequently, the bias held is reversed from to — adi- 
abatically. Then a vortex with the winding number 2F 
is formed, were F is the quantum number of the hyper- 
hne spin of the condensate. This theoretical proposal 
has been demonstrated experimentally by several groups 
IMi. A part of the condensate is lost while vortex for¬ 


mation takes place. This is due to the fact that the gap 
between the WFSS and the neutral state (NS) and the 
strong held seeking state (SFSS) is gradually closed near 
the center of the condensate during the vortex formation 
and the atoms in the trapping WFSS make transition to 
the NS and the SFSS through the Majorana hops. It 
is the purpose of this Letter to show that it is possi¬ 
ble to suppress the Majorana hops so that more atoms 
are left in the trap when a vortex is formed. It is cer¬ 
tainly desirable to have more atoms left in the trap since 
higher atomic density implies (i) easier measurement of 


physical quantities (ii) hghting against atom loss due to 
hyperhne spin changing collision and (iii) smaller vortex 
core size, which is advantageous for for ming a vortex in 
an arbitrary position in the condensate [l6|. Motivated 
by the concept of shortcuts to adiabaticity 0 , we use 
the counter-diabatic formalism (or quantum transition¬ 
less algorithm) developed by Demirplak and Rice 0 and 
formulated by Berry [Tol . 0 and demonstrated exper¬ 
imentally in [21| . The supplementary counter-diabatic 
held (CDF) thus obtained, is, however, unphysical in that 
it violates the Maxwell equation div B = 0 and some 
approximation must be introduced so that the counter- 
diabatic control is physically feasible. The approximated 
CDF turns out to be a quadrupole held, which can be 
implemented with an ordinary Ioffe bars. However, the 
CDF is generated by another set of Ioffe bars, which is 
obtained by rotating the original conhning Ioffe bars by 
7r/4. This means that two sets of lohe bars are required to 
physically realize this nonadiabatic control, which makes 
its experimental demonstration rather demanding. To 
overcome this difficulty, we apply the gauge transforma¬ 
tion [ 2 ^ to make the combined magnetic held (the orig¬ 
inal conhning held and the CDF superposed) parallel to 
the original conhning quadrupole held. The modest price 
one has to pay by this gauge transformation is to modify 
the time-dependence of the bias held B^it) from a lin¬ 
ear behavior. Experimental demonstration is made much 
easier with this modihcation. 

We hrst regard each atom with hyperspin F as an inde¬ 
pendent quantum system hxed at r and consider nona¬ 
diabatic control thereof. Let B{r,t) = {B±{r), Bz{t)) 
be the local magnetic held with the quadrupole held 
B±{r) = B'j_{x,—y) and the uniform bias held Bz{t) = 
B^{1 - 2t/T) for 0 < t < T, B^(t) = -B^ifi) for t > T 
where T is the inversion time of Bz{t). An atom at (r, t) 
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interacts with the magnetic field through the Hamilto¬ 
nian 


HBir,t) =-fB{r,t) ■ F, (1) 


where 7 = fJ-ngp, Mb is the Bohr magneton, gp is the 
g-factor for the hyperfine spin F and F = (F^, Fy, F^) 
is {2F -|- l)-dimensional irreducbile representation of the 
Lie algebra su(2). We take F = 1 in the following for 
definiteness. For a given the Hamiltonian ([T]) 

has three normalized eigenstates 
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where B = ^B'\^{r) -|- B'l{t) and (f) is the azimuthal 
angle. The corresponding eigenvalues are | 7 |i?, 0, and 
— |7|i?, respectively, where 7 = —mb/ 2 . The WFSS with 
the highest energy lyjH is the only state that can be mag¬ 
netically trapped. We start with a pure WFSS at t = 0 
and want to keep as many atoms in the WFSS as pos¬ 
sible at t = T, when the vortex formation is complete, 
in spite of vanishing gaps among the three states at the 
origin (r = 0) at t = T/2. This is realized by adding a 
CDF on the atom. 

Here we briefly summarize the counter-diabatic ap¬ 
proach to nonadiabatic quantum control and apply 
it to our hyperfine spin control. Let Ho{t) be a 
time-dependent Hamiltonian and \n(t)) be an instan¬ 
taneous eigenvector with the eigenvalue En such that 
Ho{t)\n{t)) = En{t)\n{t)). In the adiabatic approxi¬ 
mation, the solution to the time dependent Schrddinger 
equation is |' 0 n(t)) = |n(t)), where 


7n(i) = ^ [ dt'En{t')-i [ dt'{n{t')\dt'n{t')). (3) 

“ Jo Jo 

Let U{t) = J2n l^'i(0)(^(0)l be the time-evolution oper¬ 
ator such that 17(t) : |n(0)) !->■ \ipnit)). The operator [/(t) 
defines a Hamiltonian H{t) = ih{dtU{t))W[t] for an ar¬ 
bitrary time evolution, not necessarily adiabatic, namely 
\'4’n{t)) is the exact solution of 


ihdt\'il)n{t)) = H{t)\il)n{t)), (4) 


where |^„(0)) = |n(0)) is satisfied by definition. If 
we write H{t) = Ho{t) + HcDit), the counter-diabatic 
Hamiltonian EJqd (t) is written as 

(5) 

n 


when the “parallel” condition, {n{t)\dtn{t)) = 0 , is satis¬ 
fied. 

Let us apply the counter-diabatic scheme to our hy¬ 
perfine spin system by taking Hq and |n) as Ho(r,t) = 
and the instantaneous eigenstates in Eqs. ([2|). 
Note that the coordinate r here is just a parameter spec¬ 
ifying the position of the atom in the condensate, and 
the geometric phase in Eq. ([3]) vanishes. By substitut¬ 
ing Eq. m into ®, the counter-diabatic Hamiltonian is 
obtained as 


Hcd = 7 -Bcd • F, ( 6 ) 


where 

BcD(r,t) = ^ y^^ (M,x,0) (7) 

7 T B^{r,t) 

is the CDF. This magnetic field looks as if it could be 
generated by four Ioffe bars that are obtained by rotat¬ 
ing the original Ioffe bars by an angle 7 r /4 around the 
z-axis. An important remark is in order here. When the 
counter-diabatic scheme is applied to a quantum system, 
it produces a counter-diabatic potential (the magnetic 
field in our case) that prevents the quantum system from 
escaping from an adiabatic time evolution. There is no 
guarantee, however, that the counter-diabatic potential 
will be physically feasible. In fact, Bqd in Eq. 0 does 
not satisfy divS^D = 0- This is not unexpected since we 
did not consider the Maxwell equation while we obtained 
Bcy)- Tbe way out of this problem is to fix the coordi¬ 
nate r to ro in B^(r,t) in the denominator of B{r,t) so 
that Bci){,r,t) oc (y,x,0) and divScD vanishes. In our 
calculation below, we show that the CDF increases the 
number of atoms forming the vortex in a wide range of 
the inversion time T by using several values of rg. The 
performance may be further improved if we optimize the 
number of atoms left in the trap by varying rg, which 
will be worked out in our future publication. We call the 
combined magnetic field as B = B + Bcu ■ 

Now we solve the Gross-Pitaevskii equation (CPE) 
with the designed magnetic field B, 

ih3i^ t) = 4- gjidmn ^ ^ l^pl 

p 
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where 1 , m,n,p £ {— 1 , 0 , 1 }, a € {x, y, z}, 
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with the mass of the atom M, and the chemical poten¬ 
tial M (the eigenvalue of the CPE at t = 0). To begin 
with, we need to find the initial condition to solve the 
CPE ([5]) . Assuming that the initial state is in the WFSS 
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|4'(0)) = /(r)|WFSS(0)), the condensate wave function 
is obtained by solving the stationary GPE 
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where g = gn + gs, B{r,t) = {[-B_l(t') +BcD(r,t)]^ + 
Bl{t)}^/^, and the WFSS is parametrized as 
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FIG. 1. (Color online) Time dependence of the ratio 
|i?CD|/|Sn|, which is independent of r and depends only on 
t. Inset shows the bound state wave function obtained by 
solving Eq. © for ^^Na. 


Here we note that /(r) and \B±{r) + BcD(f,0)| are 
in fact the functions of r only. It turns out that Bqd 
is negligibly small at t = 0 and it can be ignored 
safely in solving Eq. ®. The time-dependent GPE is 
solved numerically and we summarize the results below. 
Throughout our calculation, we have taken the parame¬ 
ters of 23Na atoms 0, = O.OSTSa^Q&u, gs = 0 and 

Bz(0) = 1 G, B'j^ = 300 G/cm for which the harmonic 
oscillator length is found to be ano ~ 9.14 x 10“^ pm 
and Hlj ^ 3.49 x 10“^^ erg. Lengths and energies 
are scaled by these parameters. In addition, we find 
ujL{r = 0,t = 0) ~ 4.40 X 10® rad/s. The time scale 
r = 27r/a;L(0,0) ^ 1.43 ps is a reasonable measure 
of adiabaticity. The chemical potential is found to be 
p — 3.66 in dimensionless units. 

Figure 1 shows the time dependence of the ratio 
\Bcu\/\B±\ for logio(P/r) = 1.4 and ro/ano = 2, which 
is independent of r and depends only on t. The in¬ 
set is the profile of the condensate at t = 0, which is 
obtained by solving Eq. ([9]) numerically. The units of 
time, length and / are r, oho and respectively. 

Now the GP equation dS]) is solved to find interesting 
quantities. Let us show first how much atoms are left 
in the trap after the vortex formation takes place. To 
mimic the trap loss, we multiply the order parameter by 
h{r) = 0.5[1 — tanh((r — ri)/A)] with ri = dOoHO and 
A = 2aHO at each time step of the numerical simulation 
Q. In Fig. 2 we show the snapshots of the amplitude of 
the order parameter for 0 < t < lOOT with the 

GDF with the same parameter set as Fig. 1. The con¬ 
densate in the trapping region at t = lOOT is mostly in 
WFSS, while SFSS and NS are gradually removed from 
the trap in the domain r > ri due to h{r). The zero of 
the order parameter at r = 0, t = lOOT is a manifesta¬ 
tion of the vortex of the WFSS. Fig. 3 shows the fraction 
of atoms N{t)/N{0) left in the trap at t T. Since we 
replaced an unphysical Bcb with a physical Bqd satisfy¬ 
ing div Bcu = 0, the result depends on the parameter tq 



FIG. 2. (Color online) Snapshots of the amplitude of the 
order parameter for 0 < t < lOOT with the same parameter 
set as that in Fig. 1. Time intervals between snapshots are 
equal with each other. The dashed and thick solid curves 
correspond to t = 0 and t — lOOT, respectively. The unit of 
is a-l. 



FIG. 3. (Color online) Fraction N{t)/N{0) of atoms left in 
the trap long time after a vortex is formed. Since t ^ T, the 
condensate is in the pure WFSS. o shows the fraction without 
Bcu while other symbols show the fractions with Bcd for 
different ro. 
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FIG. 4. (Color online) Number of atoms N{t) for t » T in 
the control with the CDF normalized by N{t) in the control 
without the CDF. The ratio is larger for small T. 



FIG. 5. (Color online) Instantaneous fraction Nw{t)/N{0) of 
atoms in WFSS during vortex formation for logjQ(r/r) = 1.4 
and ro = 2aHO. The dashed and solid curves correspond to 
the controls without and with the CDF, respectively. 

in the denominator of Eq. O- Observe the prominent 
improvement in the ratio N{t)/N(0) for small T. This 
is more clearly seen by plotting N{t) in the control with 
the CDF normalized by N(t) in the control without the 
CDF. Figure 4 shows the ratios with the same rg as Fig. 
3. The ratio is almost 20 for T = lOr and tq = O.SaHO- 

Next, we show the CDF really suppresses the transi¬ 
tions WFSS —>■ NS and WFSS —>■ SFSS by looking at 
the projected atom numbers. For this purpose, we define 
the projection operators IIw = |WFSS)(WFSS|,IIn = 
|NS)(NS| and IIs = |SFSS)(SFSS| and evaluate the num¬ 
bers of atoms in these states as 

fVw = (d'|nw|4'),A^N = (4'|nN|d'),iVs = (^insl^-). 

Figure 5 shows the time dependence of Nw in the con¬ 
trols without CDF (dashed curve) and with CDF (solid 
curve). The figure clearly shows that the nonadiabatic 
transitions at t ^ T/2 are greatly suppressed and tran¬ 
sitions take place only at around t = T/2. This means 
that our approximation introduced in the denominator of 
(O is not justifiable at t ^ T/2, where the original Bcu 
diverges at r = 0. Otherwise, our simple approximation 


works reasonablly well. 

Now an important observation is in order. In our orig¬ 
inal proposal, we need to prepare two sets of Ioffe bars to 
produce the original confining quadrupole field and the 
counter-diabatic magnetic field. Clearly this is demand¬ 
ing for experimentalists. All topological vortex forma¬ 
tion experiments so far were conducted with a single set 
of Ioffe bars. Even if one could build a trap with two sets 
of Ioffe bars, aligning their centers exactly at the same 
place is practically impossible. To circumvent this prob¬ 
lem, we introduce time-dependent gauge transformation 
so that the combined field B± -f Ben is rotated and the 
resulting field is in parallel to the original B±. Let 

\B\ 

otB = tan“^ ——I (10) 

\B\\ 

be the angle between B and JBy. If B is rotated by —as, 
it becomes parallel to 15and hence it can be generated 
by a single set of Ioffe bars. This rotation is implemented 
by the unitary transformation 

= ( 11 ) 

Now the Zeeman term of the Hamiltonian is transformed 
as 

U{aB)lB ■ FU\aB) = "fB± ■ F, (12) 

where .Bj. = |.B|(a;, —y, 0). Of course, this is not the 
whole story and there is a price we need to pay. The 
Hamiltonian in the rotating frame acquires a gauge term 

-iUdtU^ = -aBF,. (13) 

This term, proportional to F^, works as a magnetic field 
in the 2 -direction and the bias field Bz{t) = — 

2t/T) is replaced by 

B,(t) = B,(0)(^l-|) -dB. (14) 

Note that as depends on time but not on space coordi¬ 
nates. Numerical caluculation shows that the number of 
atoms in WFSS in the control with .Bj_ and B^ is exactly 
the same as that for two sets of Ioffe bars for t >T. Fig¬ 
ure 6 shows Bz{i) and \B/\B±\ for ^^Na with the 
parameters used in Fig. [51 

In summary, we have proposed a method to suppress 
nonadiabatic transitions while topological vortex forma¬ 
tion takes place in BEC of alkali atoms. The counter- 
diabatic field is generated by a set of Ioffe bars, which is 
obtained by rotating the original Ioffe bars producing the 
confining quadrupole field by 7r/4. Our numerical calcu¬ 
lation demonstrates that nondaiabatic transtions are sup¬ 
pressed for any inverstion time T and, in particular, it is 
most impressive for a small T. We can further improve 
this scheme by applying a gauge transformation to a ro¬ 
tating frame so that the combined field B is parallel to 
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FIG. 6. (Color online) Time dependence of Bz(Q) (solid 

curve) and B^ (t) (dashed line) for the same parameters used 
in Fig. [5] The vertical and horizontal axes of the inset are 
\B±_ {r, t)\/\B± (r)| and time with the range T/A < t < 3r/4, 
respectively. 


B±. This requires modulation of Bz (t) from linear time- 
dependence. We believe our proposal is experimentally 
feasible by simple modifications of the existing setup. 
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